BCS-like Modewise Entanglement of Fermion Gaussian States 

Alonso BotercQ 

Departamento de Ftsica, Universidad de Los Andes, 
Apartado Aereo 4976, Bogota, Colombia 

Benni Rezniljt 

' Department of Physics and Astronomy, 

o" 

CN ■ Tel-Aviv University, Tel Aviv 69978, Israel. 

0) '■ (Dated: February 1, 2008) 

Q 

We show that with respect to any bipartite division of modes, pure fermion gaus- 

' 

■ sian states display the same type of structure in its entanglement of modes as that 

t^. ■ of the BCS wave function, namely, that of a tensor product of entangled two-mode 

squeezed fermion states. We show that this structure applies to a wider class of 

o 

"isotropic" mixed fermion states, for which we derive necessary and sufficient condi- 



tions for mode entanglement. 
The nature of many-body entanglement in various solid-state models such as spin chains, 



qh superconductivity, and harmonic chains, has been a recent subject of intense investigation 

[1-12]. The motivation for this effort is two-fold: on the one hand, such systems exhibit a 



rich entang 
processing 



ement structure which could potentially be harnessed for quantum information 
; on the other hand, this effort promises to deepen our understanding 
of certain universal features in many-body physics-such as quantum phase transitions- and 
their relation to entanglement [2-5,121. 

n 

In this Letter we investigate the properties of bipartite fermion mode entanglement [1( 
for a wide class of models describing correlated fermion systems, and show that this en- 
tanglement displays a simple, universal structure. As it is well known, under certain field 
transformations or after suitable approximations, a large class of interacting theories can be 
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mapped to effective theories described by quadratic Hamiltonians of the generic form 

h = CiAh + £ + h - c -) > (!) 

which are then diagonalized through Bogoliubov-Valatin (i.e., canonical) transformations}^ 
to appropriate quasi-particle mode bases. Important classes of such theories include the 
Hubbard Model and the BCS theory of superconductivity in the 'Mean Field' or Hartree- 
Fock approximation, as well as certain exactly solvable spin-chain models after a Jordan- 
Wigner transformation 11, As we shall shortly clarify, an important feature of the above 



class of Hamiltonians is that not only the ground states (the quasi-particle vacuum), but also 
all eigenstates describing a definite set of quasi-particle excitations, belong to the general 
class of fermion gaussian states. The bipartite entanglement structure of the eigenstates 
of (P) is then given by the following result, which is applicable to all pure gaussian fermion 
states, and which constitutes the main result of this Letter: 

Given a collection of N -fermion systems or "modes", partitioned into two arbitrary sets, 
A = {Ai, . . . , A m } and B = {B\ . . . , B n }, of sizes m and n = N — m, respectively, any pure 
fermion gaussian state ab of the modes may always be written as 

l^ = l^ lS J^-..|&> iS .|0>^|0>_ p (2) 

Here, s < min(m,n), and A = {Ai...,A m }, B = {Bi,...,B n } are new sets of modes 
obtained from A and B, respectively, through local fermion canonical transformations. The 
states j^fe) are two-mode fermion squeezed states of the form 

|^)_ fc - & =cos^|00)_ fe - fc -sin^|ll)_ fc _ fc) (3) 

entangling the modes A^ and Bk for k < s, and |0)~ and |0)^ are products of vacuum 
states for the remaining modes in A and B respectively. 



Together with a similar decomposition obtained for boson gaussian states in |13l . 114 1 , the 
above result shows that the bipartite entanglement structure of pure multi-mode gaussian 
states is of 1 x 1-entangled mode pairs independently of the statistics, the partition of 
modes, or the nature of the basis in which the partition is performed. The present result 
shows that in the fermion case, the mode-pairing is similar to that of the BCS wave function 
|BCS) = Uk(cos9k + sin^fcfej^&L^lvac), except that the vacuum state and the pairwise 
excitations in (j2J) are all relative to a quasi- particle spectrum determined by the multi- 
mode state |^) and the choice of partition. 



As in the boson case investigated in [13( , it can be shown that the above theorem and the 
decomposition (J2J), follows for an arbitrary pure Gaussian Fermi state from the properties 
of the Schmidt decomposition, and from the comparison of the reduced density matrices in 
their normal form, at regions A and B. In this letter we shall, however, follow a more gen- 
eral framework, which enables us to generalize this theorem to a certain class of "isotropic" 
mixed gaussian states, characterized by a symmetry property of the covariance matrices. 
The remainder of the Letter is thus structured as follows: we first review fermion gaussian 
states, canonical transformations, and define the fermion covariance matrix; next, we define 
isotropic gaussian states and prove a general decomposition for isotropic states in the lan- 
guage of covariance matrices; finally we also derive a necessary and sufficient conditions for 
the entanglement of isotropic mixed states and close with some remarks on the applicability 
of our results. 

Let a system of N fermion modes be described by a set of creation and annihilation 
operators b\,bi satisfying anticommutation relations {b\, &]•} = {bi,bj} = 0, = 
We define a fermion gaussian state for such a system as any state p that, with a certain 
choice of mode basis bi = ujbi + v^bj preserving the fermion anticommutation algebra, 
acquires the form 

N 

p = ®p k , p fc = -(i-M&jA]) , (4) 
fe=i 

with | Aj| < 1 and |Aj| = 1 for pure states (note that p can also be written in the form 
p = Z^ 1 exp(— /Sfcjfej)). Fermion gaussian states share with their boson counterparts the 
property that correlation functions for the creation/annihilation operators are completely 



determined by the two-point functions according to Wick's theorem 



16j|. Moreover, since this 



property is extensible to correlation function pertaining to a reduced subset of the modes, 



it follows that any partial (reduced) density matrix obtained from p remains gaussian |17|. 

The gaussian nature of a state is preserved under any unitary transformation that induces 
a canonical linear transformation of the fermion variables, in other words, that preserves the 
anticommutation relations. To determine the group of canonical transformations, it becomes 
convenient to replace the N creation and annihilation operators by 2N hermitian "fermion 
quadrature" combinations 

7«-i = bi + b\ , 72i = i(bi - b\) . (5) 



4 



The anticommutation relations can then be seen to be a consequence of an M. 2N Clifford (or 
Dirac) algebraic 

{la, 7/3 } = 28 a p (6) 

satisfied by the j a . Inspection of © then shows that any linear transformation of the form 
la = Oafs'jp , where O E 0(2N) (the full orthogonal group in 2N dimensions), preserves the 
Clifford algebra © and hence the anticommutation algebra of the fermion fields. 

That the group 0(2N) of fermion canonical transformations can be implemented unitarily 
on the Fock space (i.e. i' a = U'j a U), follows from the following: for any real vector v G M. 2N , 
define f — f* — v - 7, and \v \ 2 = v ■ v (= jj 1 ) with the usual inner product. We then verify 
that any real vector a with \a\ = 1 defines a unitary operator U = </l, the action of which is 
an inversion (det(O) = —1) in the orthogonal subspace to a 

flint = [Za^dy - b^Ytv ; (7) 

in IR 2 " - , this inversion is equivalent to an SO(2N) rotation times a reflection along a. Since 
any orthogonal transformation can be generated from a product of a certain number of such 
reflections along different axesj^j], it follows that any orthogonal transformation of the 7's 
can be induced by a suitably chosen unitary operator of the form U = ^1^2(^3 • • • tp with 
\a k \ = 1. 

The group of operations thus defined is the Pin(2N) group, the double cover of 



0(2N) 2l|, |22|, |23|. In turn, the subgroup of Pin(2N) generated by an even number 



of reflections is the covering group of SO(2N), Spin(2N), associated with the linear 



or squeeze 



fermion canonical transformations (also known as the Bogoliubov-Valatin^ 
transformations 26]). More generally, Pin(2N) includes improper orthogonal transforma- 
tions (det(O) = — 1) generated by an odd number of reflections, and which may in fact 
be forbidden due to fermion number superselection rules in a fundamental fermion theory. 
This fact, however, does not affect the generality of the results presented here. Since an 
improper transformation is a SO(2N) rotation times some canonical reflection, an improper 
transformation is equivalent to a proper transformation but with the role of the creation 
and annihilation operators exchanged for at most one mode. 

An important consequence of the 0(2N) equivalence of fermion gaussian states is that 
all multiparticle eigenstates obtained from a given vacuum are gaussian. Indeed, if |vac) 
is the ground state of a certain Hamiltonian, annihilated by the 6« operators in a given 
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quasi-particle basis, then for each mode it follows that b\\va.c) = 72«|vac); it therefore follows 
that multiparticle states are obtained from some U G Pin(2N) acting on |vac). Since the 
vacuum state is gaussian, the resulting state will also be gaussian. 

We next define for any gaussian quantum state the fermion Covariance Matrix (FCM) as 

^=^<[7«,7/j]>„- (8) 

Being an antisymmetric 2N x 2N matrix, the FCM can always be brought to the block 
diagonal form 

N f -1 \ 

W = OMO T = 0A,J 2 , J 2 =| I, (9) 

by means of a 0(2N) transformation, and with the Aj > 0. In terms of the creation/ 
annihilation operators obtained from the transformed gamma's 7 = O7, the gaussian state 
p then takes the form (0J with W as its FCM. By setting all the \ positive or zero, we 
obtain a unique characterization of a fermion gaussian state (up to relabellin g of the modes) 
analogous to the so-called Williamson normal form of the boson case lH Il8^. We shall 
therefore term the canonical form W as the fermion Williamson form, and the {\{\ (A, > 0) 
in W as the Wiliamson eigenvalues. 

Now, the Williamson spectrum can be obtained from the doubly degenerate spectrum of 
the matrix — M 2 . This provides a simple characterization of pure gaussian states in terms 
of their FCM; since such states are unitarily equivalent to pure states of the form with 
Aj = 1 , Vz, the pure FCM s satisfy the condition — M 2 = 1. Pure gaussian states can then 
be viewed as the extremal states in a more general class of isotropic gaussian states, for 
which the square of their FCM's remains invariant under any 0(2N) transformation: 

M 2 = -X 2 t, A <1. (10) 

It is then possible to show that if M is the FCM of an isotropic gaussian state, and if 
the modes are separated into two sectors A and B (7 = 7^ © 7^), then there exist local 
orthogonal transformations •ja = OaIa and 7b = ObIb such that the upon re-ordering the 
modes the FCM takes the form 

M = M_ _ © M_ _ © ... @M_ _ © , (11) 

A 1 B 1 A 2 B 2 A S B S A F Bp 



6 



7~ and of the form 



M 



k + A 



\ 2 



(12) 



where _ is isotropic for the entangled pair sector 7. 

/ -A, \ 

\ Ki 

— —A,, 
Xi / 

and and M are isotropic FCM's with Williamson eigenvalue Ao for the remaining 

A F Bp ' 

modes in the A and B sectors respectively. From the correspondence between FCMs and 
gaussian mixed states, it then follows that in the new mode basis, the isotropic fermion 
gaussian state p^ B takes the form 



P 



(0) 



P 



A 1 B 1 



P: 



(13) 



where p~ ~ are two-mode fermion squeezed states with FCM of the form ()12)) . and 

pio) 

Ap 

and pi> are local isotropic states for the remaining modes in A and B. We note that for an 

Bp 

FCM of the form ()12j) . the corresponding two-modes squeezed state may be written as 

l-AoPLX]W 



\ Ti ( I 



p AlB , '■ ' V' 1 . A 

with Tj an entangling unitary operator 



T; 



cxp 



— f fotfot +6-6- 



tan 29a 



K 1 



(14) 



(15) 



Thus, in the pure case Ao = 1, the two-mode entangled states become projection operators 
onto two-mode squeezed states Tj|00)- _ which expand out to the form ©, whereas the 

A i B i 

unentangled states become projection operators onto the vacuum states of the remaining 
modes. 

To prove the decomposition (|11|) . first perform local orthogonal transformations 7a©7b — 
{0 A © Ob)ia © 1b bringing each of the local FCMs M A = ±([ 7x , 7J]) , M B = ±{[ib , 
into the canonical Williamson form ©, with all \ > 0. The total FCM thus obtained may 
be written as 



M 



(16) 



W A K 
-K T W B 

with Wa = ©iliA~ J2 and Wb = ©™=iA~ J 2- Substituting into the definition of the 
isotropic matrix M 2 = — Aq, the following equation is obtained from (|16j) 



W A K + KW B = 



(17) 
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Consider then a 2x2 sub-block = — 2(7- 7 T ) of K connecting the modes with eigenvalues 
A- and A. . From equation (fTTj) and using J| = — 1 2 , we find that 

A,- 



A-^ = A- J 2 ^J 2 . (18) 

For A- 7^ A^ this equation has no solution other than Kij = 0, meaning that modes in A 
and modes B with different Williamson eigenvalues decouple. 

Thus, let 7- and 7- stand for the modes in A and B with the same local Williamson 
eigenvalue A, and group the modes according to their eigenvalues so that M takes the 
Jordan form M = @ A M x where each block M x is the FCM for the modes in A and B with 
a common local Williamson eigenvalue A. Concentrating on a given A, let 9a and 9b be the 
degeneracies of A in the Williamson spectra of Wa and Wb respectively, so that M\ may be 
written as 

\ a aj 2ss y i=i 

Now define ka = -\/Ag — A 2 and note that Ma is also an isotropic FCM with symplectic 
eigenvalue A . Thus, M 2 = — Aq yields 

K x Kl = n\t 2gA , K\K\ = K\t 2gB (20a) 
J2 9A K x J 2gB = K x , (20b) 

where we have used the fact that J 2g = —$-2g- Taking the trace of both equations in ()20a|) 
and using Tr^^A-^A ] = T^bIK^Kx], we find that K 2 x (gA — 9b) = 0, proving that for A 7^ A , 
the degeneracies of A in the local FCM's are the same. Thus, let qa = 9b = 9 and express 
K in terms of an unspecified matrix Q\ according to 

N f l\ 

K x = kQ x (3, /3 = I 1 . (21) 
Noting that [3 2 = l 2g and J 2g } = 0, we find from Eqs. (J2U|) that Q\ satisfies 

Q\Q\ — Q\Q\ — l2 5 QxJ2gQ\ 1 — J2g , (22) 

indicating that Q\ is an orthogonal symplectic transformation. Hence, K can be brought 
to the block diagonal form K' x = K\f3 by means of a one-sided orthogonal transformation 
O a = <5l' leaving the upper diagonal block of (|T§j) intact. Grouping the resulting modes 
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in each degenerate space as 7a = ©f =1 (7 A ©7~ ), we achieve the pair- wise decomposition 
for A 7^ Ao, with each pair described by a covariance matrix of the form (jl2j) . For A = Ao, 
the degeneracies on each side are not restricted, and equations (ffilj) imply that K\ Kj Q = 
=>- K\ = 0. Therefore, local modes with symplectic eigenvalue Ao decouple, as may be 
expected for the pure case Ao = 1. 

We briefly comment on the quantification of bipartite entanglement as obtained from the 
modewsie decomposition. In a fundamental fermion theory, the characterization of entan- 
glement is somewhat delicate, as one needs to take into account the absence of a natural 
tensor-product decomposition of the fermionic Hilbert space, as well as restrictions on the 
set of possible local operations posed by the anticommutativity of fermion operators |22t and 
superselection rules due to fermion number conservation 28]. Two measures are available: 
the explicit entanglement content of the state with respect to the second quantized mode 
Fock space-quantified by the "entanglememt of modes" E M and the operationally ac- 
cesible bipartite entanglemement of the state-quantified by the "entanglement of particles" 
Ep, with Ep < Em [28]. For the pure case of Eq. 0, Em will simply be the sum of the von 
Neumann entropies of the partial density matrices obtained from the Schmidt coefficients 
of the individual entangled mode pairs in Eq. ©• On the other hand, a characterization in 
terms of Ep is less direct, as the measure is superadditive and will depend on the funda- 
mental fermion content of the mode bases in the decomposition. Nevertheless, since Ep and 
Em coincide asymptotically j3], one may expect that for entangled states involving a large 
number of fermions, Em should give a reasonable indication of the usable entanglement. 

If mixed fermion states are analyzed according to the Em criterion, it is possible to asses 
the separability of isotropic gaussian states via the Peres-Horodecki partial transpose (PT) 
criterion ^| . Since a two-mode mixed state (pPfj) is equivalent to a two-qubit mixed state, 
PT-negativity of at least one of the states p_ - in Eq. (H3J) will then be sufficient £0j to 

A s B s 

ascertain the inseparability of the isotropic state p(°\ Using the fact that the entangling 
operator Tj in (fTK|) performs a rotation between |00)- - and |11)~ - while leaving |01)~ - 
and 1 10) ~ - unaltered, the matrix representation of pi can easily be shown to take the block 



form: 




m^\ MX ' , +Ki 2K : L r~ A ° „ i P3) 



2«j (1-A,) 2 + ^ J \ 1-Ag 
where the first block refers to the 1 00) - - and |11) - ~ vectors and the second to the 1 01) 
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and 1 10)- ~ vectors. In this decomposition, the partial transpose operation amounts to 
swapping the off-diagonal elements between the two blocks. The negative partial transpose 
criterion for pi then becomes Ao > ^% > |(1 — Aq). Note that this condition cannot be 
satisfied for any Kj if Ao < V% — 1 and is satisfied for all non-zero K; t in the pure case Ao = 1 
as expected. 

In conclusion, we have studied fermion mode entanglement and showed that under an 
arbitrary bi-partite division of a set of fermion modes, gaussian mode entanglement has a 
universal structure of a generalized BCS-like wave-function. Interestingly, unlike the gaus- 
sian boson case, our result applies not only to the ground state, but also to all definite 
excitation eigenstates of quadratic fermion Hamiltonians. Our result can be used for com- 
puting the bi-partite mode entanglement of a variety of states in a wide class of models, 
such as certain spin-chain models, and fermion liquids and superconductivity in the mean- 
field approximation, we believe that the present results can be helpful in connection with 
information processing applications involving coupled many-particle fermion systems or su- 
perconductors. For instance, when two such systems are coupled by a quadratic tunnelling 
interaction [3] , our result implies that under a proper choice of "local operators" , the total 
wave function factorizes to a product of non-maximal EPR-like states. 

A.B. acknowledges support from Colciencias (contract 1204-05-13622). B.R. acknowl- 
edges the support of ISF grant 62/01-1. 



[1] P. Zanardi, Phys. Rev. A. 65 042201 (2002). P. Zanardi, X. Wang J. Phys. A, 35 7947 (2002). 

[2] K Audenaert, J. Eisert M. B. Plenio and R. F. Werner, Phys. Rev. A 66, 042327 (2002). 

[3] M.B. Plenio, J. Hartley, J. Eisert, New J. Phys. 6, 36 (2004). 

[4] V. Vedral, |quant-ph/0302040| 

[5] F. Verstraete, M. Popp, J.I. Cirac, Phys. Rev. Lett. 92, 027901 (2004). 



[6] A. Botero and B. Reznik, quant-ph/0403233 



[7] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120126 (1998). 

[8] P. Recher and D. Loss, Phys. Rev. Lett. 91, 267003 (2003). 

[9] C.W.J. Beenakker, M. Kindermann, Phys.Rev.Lett. 92, 056801 (2004) 

[10] F. Verstraete and J.I. Cirac, |quant- ph/0311130 



10 



[11] J. Gonzalez, M.A. Martin-Delgado, G. Sierra and A. H. Vozmediano, "Quantum Electron 

Liquids and High-T c Superconductivity", (Springer, 1995). 
[12] E. S. Fradkin, "Field theories of condensed matter systems" (Addison- Wesley, 1995). 
[13] A. Botero and B. Reznik, Phys. Rev. A 67, 052311(2003). 
[14] G. Giedke, J. Eisert, J.I. Cirac, M.B. Plenio Quant. Inf. Comp. 3, 211 (2003) 
[15] M. C. Chung and I. Peschel, Phys. Rev. B 64, 064412 (2001). 
[16] M. Gaudin, Nucl. Phys. 15 89 (1960). 
[17] I. Peschel, J.Phys.A: Math.Gen. 36, L205 (2003) 

[18] J. Williamson, Amer. J. Math. 58, 141 (1936); R. Simon, E.C.G. Sudarshan and N. Mukunda, 

Phys. Rev. A 36 3868 (1987). 
[19] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[20] M. Horodecki, P. Horodecki and R. Horodecki, Phys. Lett. A 223, 1 (1996). 
[21] P. Lounesto. "Clifford Algebras and Spinors", (Cambridge, 1997) 

[22] Y. Choquet-Bruhat and C. DeWitt-Morette, "Analysis, Manifolds and Physics" , Vol. 2. (North 

Holland, Amsterdam, 1989) 
[23] I. M. Benn and R. W. Tucker. "An Introduction to Spinors and Geometry", (Adam Hilger, 

Bristol, 1987) 

[24] F. A. Berezin, "The Method of Second Quantization" (Academic Press, New York, 1966) 
[25] N. N. Bogoliubov, Nuovo cimento 7, 794 (1958); J. G. Valatin, Nuovo cimento 7, 843 (1958). 
[26] K. Svozil, Phys. Rev. Lett. 65, 3341 (1990). 
[27] S. Bravyi, A. KitaevJqua^ph/0003137 

[28] H. M. Wiseman, J. A. Vaccaro, Phys. Rev. Lett. 91, 097902 (2003). 

[29] M. H. Cohen, L. M. Falicov and J. C. Fillips, Phys. Rev. Lett. 8, 316 (1962). 



